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We investigate the evolution of field line helicity for magnetic fields that connect two boundaries without 
null points, with emphasis on localized finite-B magnetic reconnection. Total (relative) magnetic helicity is 
already recognized as an important topological constraint on magnetohydro dynamic processes. Field line 
helicity offers further advantages because it preserves all topological information and can distinguish between 
different magnetic fields with the same total helicity. Magnetic reconnection changes field connectivity and 
field line helicity reflects these changes; the goal of this paper is to characterize that evolution. We start 
by deriving the evolution equation for field line helicity and examining its terms, also obtaining a simplified 
form for cases where dynamics are localized within the domain. The main result, which we support using 
kinematic examples, is that during localized reconnection in a complex magnetic field, the evolution of field 
line helicity is dominated by a work-like term that is evaluated at the field line endpoints, namely the scalar 
product of the generalized field line velocity and the vector potential. Furthermore, the flux integral of this 
term over certain areas is very small compared to the integral of the unsigned quantity, which indicates that 
changes of field line helicity happen in a well-organized pairwise manner. It follows that reconnection is very 
efficient at redistributing helicity in complex magnetic fields despite having little effect on the total helicity. 


I. INTRODUCTION 

Magnetic helicity is a valuable concept in magnetohy¬ 
drodynamics (MHD) that quantifies the linking, twist¬ 
ing and kinking of magnetic field finest. It is highly 
conserved under a broad range of circumstances and it 
therefore has many applications in both laboratory and 
astrophysical plasmas^. More precisely, total magnetic 
helicity (the volume integral of A • B where A is the vector 
potential and B = V x A is the magnetic field) is an ideal 
MHD invariant for magnetically closed domains^, which 
is readily extended to magnetically open domains either 
as total relative magnetic helicity^^ or by an appropri¬ 
ate restriction of the gauge of the vector potential^. Ideal 
invariance holds because ideal evolutions neither create 
nor destroy magnetic helicity and are unable to transport 
helicity across the magnetic field. In non-ideal MHD, ex¬ 
act conservation of magnetic helicity breaks down but 
magnetic reconnection at high magnetic Reynolds num¬ 
ber nonetheless conserves total (relative) magnetic he- 
licity very well^. Reconnection does, however, redis¬ 
tribute helicity between field lines as magnetic connectiv¬ 
ities change. Thus, magnetic reconnection approximately 
conserves total magnetic helicity but may radically alter 
how that total is composed. 

An extreme example of the redistribution of magnetic 
helicity is Taylor relaxation. Considering a reversed-field 
pinch, Taylor— hypothesized that turbulent magnetic 
reconnection allows an initial magnetic field to relax to 
the minimum energy state with the same total helicity, 
which had previously been shown by Woltjer— to be a lin¬ 
ear force-free field. This assumes that total helicity is the 
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only helicity constraint and requires complete redistribu¬ 
tion of magnetic helicity across the cross-section of the 
device. Taylor relaxation was successful for reversed-field 
pinches and has since been investigated for other situa¬ 
tions e.g. for the solar corona with applications to coro¬ 
nal heating and microflaresii - — . There are also known 
examples where the end state is a nonlinear force-free 
fieldii - — , however redistribution of helicity, although less 
extensive, is a major feature of those cases as well. 

Due to its conservation, magnetic helicity is of broad 
astrophysical interest, especially in scenarios involving 
magnetic reconnection. For instance, the generation of 
magnetic fields by dynamo action is intrinsically related 
to the properties of magnetic helicityi£. To give a few 
more examples, helicity conservation has been invoked in 
magnetospheric physics to explain generation of twisted 
flux tubes during dayside reconnection and plasmoid for¬ 
mation in the magnetotai l 18 i 19 . In solar physics, mag¬ 
netic helicity is injected into the corona by flux emer¬ 
gence and photospheric motions including differential ro¬ 
tation and shearing flows in active regions^. Values and 
changes of magnetic helicity in active regions have been 
linked to solar flares and coronal mass ejections^! - — , he- 
licity “condensation” is a candidate explanation for the 
formation of filament channels^, expulsion of helicity 
from corona leads to the presence of twisted flux ropes in 
the heliosphere^ and reconnection of flux tubes can be 
a source of torsional Alfven wave a 26 i 27 . 

In this paper we consider a refined measure of helic¬ 
ity: a helicity density which is assigned to each field line. 
This “field line helicity” contains all available topological 
information and can therefore distinguish between mag¬ 
netic fields with the same total helicity. The primary aim 
is to investigate how this measure of helicity evolves in 
the broad regime between ideal evolution (for which ev¬ 
ery field line has its own helicity invariant) and Taylor 
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relaxation (for which the only helicity invariant is total 
helicity). The results characterize changes to the compo¬ 
sition of magnetic helicity and are expected to advance 
our understanding of 3D magnetic reconnection across 
a broad variety of applications including turbulent mag¬ 
netic relaxation. 

This paper is organized as follows. Section ITT1 describes 
the model, recaps the concept of field line helicity and 
discusses gauge considerations. In Sec. Hill the evolution 
equation for field line helicity is derived and its terms are 
examined. We then focus on cases where dynamics are lo¬ 
calized within the domain (Sec. El and show that evolu¬ 
tion of field line helicity for a given field line is dominated 
by a work-like term which has a well-organized structure 
of pairs of positive and negative rates of change. Kine¬ 
matic examples in Sec. [V| confirm our analytic results. 
The paper ends with a summary in Sec. |VlJ 


II. PRELIMINARIES 
A. Model 

Much of what we discuss in this paper applies gen¬ 
erally, but for concreteness we will consider a flux tube 
model sketched in Fig. [TJ All field lines enter through 
a single surface Dq and exit through a different surface 
D\. The remainder of the boundary is a magnetic sur¬ 
face, Ds , that joins the edge of Dq with the edge of D\. 
There are no magnetic null points in the domain, hence 
we consider finite-B reconnection^. This model can de¬ 
scribe closed flux tubes (in which case magnetic field is 
periodic on Dq and D\) as well as open flux tubes. More 
general magnetic fields can be partitioned into a collec¬ 
tion of such domains, which increases the generality of 
this model. 

For simplicity, some of our results will be presented 
using a restricted version of the model. The first simpli¬ 
fication is to take Do and D\ planar with outward sur¬ 
face normals, n, pointing in the indirection. In these 
cases, we also assume that the boundary is line-tied and 
ideal so that electric field E = 0 and B • n is constant 
in time. This restriction is appropriate to solar physics 
where the restricted model may, for example, represent 
the magnetic field in a coronal loop under assumptions 
that the coronal magnetic field is evolving more rapidly 
than the photospheric convection timescale and that dy¬ 
namics are concentrated away from the side boundary, 
Ds- The simplified model is therefore of direct relevance 
to the relaxation of magnetic braids, magnetic instabili¬ 
ties, coronal heating and confined solar flares. 



FIG. 1. Sketch of the flux tube model with all field lines con¬ 
necting Do to D i. The side of the domain, Ds , is a magnetic 
surface and the domain contains no magnetic nulls. 


flux with the field line of interest^ and can also be viewed 
as a topological flux function^ - —. It is defined as 

A(x) = f A • dl (1) 

J F(x) 

where x is a point on a cross-section of our flux tube 
(e.g. the lower boundary, Do), F(x) is the magnetic field 
line through that point and dl is the line element along 
the field line. Since our domain is free of null points the 
integral is always well defined. 

Field-line helicity retains all topological information, 
in contrast to the volume-integrated total magnetic 
helicity^. It can therefore distinguish between topologi¬ 
cally different magnetic fields with the same total helicity, 
while the total magnetic helicity is easily recovered from 
A since 

H(A) = [ A-B d 3 x= f ABid 2 x= f AB a d 2 x 
Jv j Do J D 1 

( 2 ) 

where Bi is the magnetic field component parallel to the 
inward normal on Do and B 0 is the magnetic field com¬ 
ponent parallel to the outward normal on D\. Note that 
this formula justifies the name field line helicity since to¬ 
tal helicity is the flux integral of A. Provided the gauge is 
suitably restricted, A has the desirable property of being 
an ideal invariant, save for changes of field line connec¬ 
tivity caused by motions on the boundaries. The gauge 
condition under which this is true is discussed in Sec. Ill Cl 
and the result is derived in Sec. IIIIBI At the same time, 
A is considerably easier to work with than the helicity 
density, A • D, which depends on all three coordinates 
and changes under ideal evolutions. 


C. Gauge considerations 


B. Field line helicity 

Field line helicity, A, assigns a helicity value to every 
field line. Physically, it measures the winding of magnetic 


Magnetic helicity and field line helicity are in general 
gauge-dependent for open domains, i.e. they change un¬ 
der gauge transformations of the vector potential. Re¬ 
ferring to Eq. a gauge transformation A! = A + Vy 
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implies a new field line helicity 

-4'(z) = J (-4 + V X )'df = A(x) + [x]^ , (3) 

where oco and x\ represent the start and end points of the 
field line respectively. It follows that field line helicity, 
although generally gauge-dependent, is invariant for a 
restricted set of gauge transformations that have x(xo) = 
x(^i). Therefore, A may be made gauge invariant by 
imposing a boundary condition on the vector potential 
that fixes the components of A tangent to the boundary. 

The physical interpretation of A and the fixing of n x A 
on dV can be regarded in two complementary ways. The 
first approach is to consider relative magnetic helicity^^, 

H R {B\B Tei ) = J d 3 ^, (4) 

where H ref = V x A ref is a reference field that satisfies 
B rei • u\qv = B • n\dv • If is common practice to choose 
the potential field for B rei . Relative helicity has the ad¬ 
vantage of being gauge invariant but has its own disad¬ 
vantage of depending on the chosen reference field, which 
may need to change in time to match B(t) on dV . If the 
gauge of A is restricted such that 

A x n\dv = A rei x n\dv, (5) 

which is always possible given B rei • n\dv = B • n|av, 
then A becomes gauge invariant as noted above. 
It can also be shown that Eq. ® and © imply 
that H(B) = H R (B\B rei ) + H(B re f). Therefore, if the 
boundary condition given by Eq. (|5j) is imposed using a 
reference field for which H(B re f) = 0, then H = Hr 
and A can be regarded as the density per unit flux of the 
gauge-independent relative helicity^i. 

An alternative approach arises from work by Prior and 
Yeates— who considered the physical interpretation of he¬ 
licity in open domains as measuring the winding between 
pairs of field lines. This can be viewed as a generalization 
of the linking number interpretation of helicity in closed 
domains^. Prior and Yeates— concluded that each gauge 
measures winding with respect to a particular frame, but 
some gauges include a non-physical contribution equiva¬ 
lent to measuring winding in a twisted frame. It is there¬ 
fore desirable to restrict oneself to gauges that measure 
winding with respect to a fixed basis, which makes A 
gauge invariant for a given frame. Furthermore, it has 
been shown that the choice of basis used to measure wind¬ 
ing of field lines corresponds directly to a choice of refer¬ 
ence field for relative helicity^, so there is an equivalence 
between the relative helicity and winding interpretations. 

The general equations derived in Sec. Hill and [IV] are 
valid in any gauge of A. In the examples presented in 
Sec. [VJ we consider situations where the reference field 
is fixed in time (which is possible since B • n is time- 
invariant on dV). Furthermore, we choose a gauge such 


that H = Hr with a potential reference field for Hr. 
Since these examples have uniform B-n on the boundary, 
this potential field is untwisted in the sense of Prior and 
Yeates—, meaning that A may be interpreted either as the 
average winding of field lines with respect to an untwisted 
basis, or as a density for the relative helicity. Since we 
have fixed the gauge of A on the boundaries over time, 
any change of A in time must correspond to a non-ideal 
process. 

III. EVOLUTION EQUATION 
A. Derivation 

The evolution equation for field line helicity can be 
derived as follows. We start from a general form of Ohm’s 
law, expressed as 

E + v x B = R, (6) 

where R is the non-ideal contribution to the electric field. 
For any R, the non-ideal term can be decomposed into 
a gradient term that produces the parallel electric field, 
plus a term perpendicular to B that sums with the gradi¬ 
ent term to give the correct perpendicular electric field, 
i.e. 

R = VT - u x B. (7) 

In this decomposition, T and u may not have simple 
analytic forms derivable from the local value of R, in fact 
they must generally be defined in a non-local manner. 
Noting that Eq. (|6j) and Eq. ([7]) imply • VT = £j| t 
is obtained from the integral 

T(f) = [ E-dl. (8) 

J F(x) 

Here, we consider situations without null points or closed 
field lines, hence T is well-defined and single-valued 
throughout the domain. Having obtained T by integra¬ 
tion of the parallel electric field, the components of u 
perpendicular to B are readily computed from Eq. 
completing the decomposition. 

Evolution of the vector potential is governed by 

f)A 

_ = (9) 

where <f> is the electric potential. Using Eq. © and © 
to substitute in (|9j), 

dA 

— = w x B — V(T + $), (10) 

where 

w = v + u 


(ii) 
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is the generalized field line velocity. 

The concept of generalized field line velocity has been 
discussed in detail by Hornig and Priest— and Hornig—. 
Briefly, taking the curl of Eq. m shows that the mag¬ 
netic field evolves as if advected with w. The slip velocity 
u = w — v (and hence w) generally depends on where one 
sets T = 0. This nonuniqueness represents the fact that 
when magnetic connectivity is changing, the generalized 
velocity of a field line depends on what plasma element it 
is traced from. Another point worth noting is that since 
u enters Eq. © as u x B, the component e# • u is not 
constrained so one may choose a single component of u 
or equivalently w. 

We are now ready to derive an equation for dA/dt. 
Starting from the time derivative of Eq. ©, we differ¬ 
entiate under the integral sign using the Lie-derivative 
formula for line integrals through 3D spaced - —, then 
simplify using Eq. m to obtain the line integral of a 
gradient, which is readily evaluated. Thus, 


dA(x,t) _ d r 

d t dt J F 

-l 

= I v( 

JF(xd) V 


A • dl 


F(x,t ) 

F(x,t) 


dA 

— - fiJxVxA + V 

at 


(^w • A^j 


w • A - T - <f> ■ dl 


w • A - T - $ 


■dl 


( 12 ) 


where xo and x\ represent the start and end points of 
the field line respectively. 

Thus, the evolution of field line helicity depends on 
the motion of field line end points parallel to the vector 
potential on the boundaries, the integral of E\\ along the 
field line of interest (i.e. the voltage drop) and the dif¬ 
ference in electric scalar potential between the field line’s 
end points. 


line helicity for relative helicity with a time-independent 
reference field^i. 

The T terms correspond to a net voltage drop along the 
field line. This quantity, f F E • dZ, has played prominent 
role in general magnetic reconnection. Such a voltage 
drop across a localized non-ideal region is necessary and 
sufficient for the change in connectivities to be felt out¬ 
side that region, i.e. a net voltage drop across a localized 
non-ideal region distinguishes finite-B reconnection with 
global effects from finite-B reconnection with only local 
effects^. Voltage drops are also widely used to measure 
reconnection, with the maximum (unsigned) voltage drop 
quantifying the reconnection rate^r— . 

Lastly, the w • A terms represent motion of field line 
end points on the boundaries. Their form is analogous to 
work done against a force and depends only on the com¬ 
ponent of w parallel to the vector potential. This term 
can be present for ideal evolutions when motions on the 
boundaries change field line connectivity, in which case 
w • A = v • A. It is also present during localized recon¬ 
nection when w is due to field line slipping. The terms 
are obviously gauge dependent since a change of A will 
in general change the pattern of w • A on Do and D i 
but they can be made gauge independent by fixing Axn 
on dV (as recommended to make the <f> terms cancel for 
every field line) and using freedom of the parallel com¬ 
ponent of u to set w • ft = 0 on Do and D\. It will also 
be shown in Sec. E22 that the contribution integrated 
over certain areas is independent of gauge even when no 
boundary condition is imposed on A. Removing this term 
entirely by gauge choice would require a time-dependent 
gauge on the boundary, which would remove the physi¬ 
cal interpretation of A as measuring changes in field line 
connectivity within the domain. 


IV. LOCALIZED RECONNECTION 


B. Interpretation of terms 

The 4> terms in Eq. m correspond to a difference of 
scalar potential between the endpoints of the field line 
and represent helicity flux along the field line due to the 
gauge. It is the only term that remains in an ideal evo¬ 
lution with no motions on the boundaries. It is usually 
convenient to use a gauge in which this term vanishes 
so that helicity flux is due to physical terms only and A 
becomes an ideal invariant save for changes of field line 
connectivity caused by motions on the boundaries. As 
an example, in our restricted model with E = 0 on dV 
one can impose a condition that Axn on dV is con¬ 
stant in time, which makes <f> spatially constant on dV 
by Eq. ©, hence the <f> terms cancel for every field line. 
Physically, this restriction ensures that A measures wind¬ 
ing with respect to a time-independent frame, eliminating 
non-physical changes^, or equivalently that A is a field 


A. Simplifications 

Several simplifications can be made when reconnection 
dynamics are localized within the domain. For simplicity, 
we do this here for planar and horizontal boundaries Dq 
and Di. We also assume that the electric field vanishes 
on the domain boundary since dynamics are internal. 

First, it is convenient to use the freedom available in 
the definition of u to set w-e z = 0 throughout the domain, 
which closes Do and D\ to transport of magnetic flux 
with w. Under this choice, 


e z x (VT - E 



which may be confirmed by using Eq. © to substitute 
for R in Eq. (|6]), then taking the cross product with e z 
and rearranging for w. Since E is assumed zero on dV, 
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the field line velocity on the boundaries is simply 


w = 


e z x W 

B z 


(14) 


which implies that field line end points move along con¬ 
tours of T. 

Next, we can choose to set T = 0 everywhere on Do, 
which gives w = 0 there by Eq. ed- We can also use a 
gauge condition that n x A is time-independent to make 
the <f> terms vanish (Sec. IIII Bj) . thereby ensuring that 
A is constant in the absence of reconnection. Doing so, 
Eq. (fl2|) simplifies to 


pi y I _^ 

— = w(xi) ■ A(xi) - T(£i), (15) 


including the vector potential^. This is because, for a 
sufficiently complex magnetic field, a pair of field lines 
traced from nearby starting points may take very differ¬ 
ent paths through the domain, thus acquiring very dif¬ 
ferent contributions to the integrand. For example, when 
the parallel electric field is integrated along field lines to 
obtain T, one field line may pass through a non-ideal 
region that the other field line bypasses altogether. In 
this way l becomes much less than L. We therefore con¬ 
clude from the scaling analysis that changes to a field 
line’s value of A occur primarily through the work-like 
term when magnetic reconnection occurs in a complex 
magnetic field. 


C. Paired increases and decreases 


which contains only terms that are evaluated at a single 
boundary (the one where T ^ 0). 

Finally, using Eq. (ED and noting that V x e z = 0, on 
Do and D\, 


w ■ A = — 


d'Vx ('he*) 




V- 


(I x (¥e*)) 


B z 

V • (’JA x e z ^j 
~Bl 


^e z - Vxi 
~B Z 


- T. 


(16) 


i.e. the work-like term can also be expressed as a diver¬ 
gence term minus the voltage drop along the field line. 


B. Dominance of work term 

The character of the evolution of field line helicity de¬ 
pends on which term dominates the right hand side of 
Eq. (USD. Referring to Eq. ED, the work-like term scales 
as 


w-A~ yT (17) 

where L is the length scale associated with the vector 
potential on the boundary ( B z = e z • V x A A/L) 
and l is the length scale associated with T on the bound¬ 
ary (||VT|| ^ T//). For the purpose of this scaling ar¬ 
gument, A and T are characteristic values. The nature 
of the evolution equation therefore depends on the ratio 
of length scales L/l, with the work-like term dominating 
for L/l 1. 

The property L/l 1 is characteristic of magnetic 
fields with complex field line mappings, for example mag¬ 
netic braids. In these cases, the complexity of the field 
line mapping means that field line integrated quantities 
such as T vary on a length scale much smaller than the 
typical gradient length scale of quantities on the bound¬ 
ary that have not been integrated through the domain, 


Having considered the change of field line helicity for a 
chosen field line, it is also instructive to look at changes 
integrated over area. Integrating Eq. (fl6j) over an area 
S C D\ and using the weight B z for consistency with the 
total helicity integral defined in Eq. (|2]), one finds 

/ B z w • A d 2 x = / TA x e z • rh dx — / B z d/ d 2 x 

Js Jds Js 

(18) 


where we have used the divergence theorem and rh is an 
outward edge normal on dS. 

There are certain choices of S for which the first term 
on the right hand side of Eq. ED vanishes to leave 


f B z w • A d 2 x = — f B z d/ d 2 x. (19) 

Js Js 

The most obvious of these is when S = D i, since T = 0 
everywhere on dD\ from the assumption that the side 
boundary of V is ideal. Thus, although the work term 
dominates changes of field line helicity for individual field 
lines, its net effect over the entire domain is the same as 
the integrated effect of the typically much smaller volt¬ 
age drop term. We conclude that the divergence part of 
the w • A term occurs as pairs of opposite polarity, which 
redistribute magnetic helicity and are individually typi¬ 
cally much stronger than T but which cancel one another 
in the area integral. 

Area integrals of the voltage drop are also readily in¬ 
terpreted. Taking the time derivative of Eq. (|2j) and as¬ 
suming that D i and B z \d 1 are fixed in time, the rate of 
change of total helicity is 


d H 
d t 



B Z A d 2 x 


L 


B,— cl c 


( 20 ) 


Then, expanding d^4/d£ with Eq. (fl5l) and using Eq. (fl9l) 
to replace the integral of B z w • A, one finds 


d H 
dt 



*) 


d 2 x = -2 



D-.T d 2 x, 

( 21 ) 
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FIG. 2. Construction of an area inside which positive and 
negative polarities of dA/dt cancel, (a) shows the curves used 
to construct the edge of the area: T = 0 contours (black) and 
lines perpendicular to the vector potential on the boundary 
(red), (b) shows an enlargement of the shaded area labeled 
to indicate which quantity is zero on each edge. It is readily 
seen that the line integral of TA x e z • fa around this shape 
is zero. 


which is equivalent to the volume integral of —2 E • B. In 
other words, —2T (half from the explicit term in Eq. (fl5l) 
and the other half implicit in the work-like term) repre¬ 
sents the net imbalance of helicity sinks and sources along 
a field line, which contributes to changing the total he¬ 
licity and is distinct from the redistribution of helicity 
caused by the divergence part of the w • A term. The con¬ 
dition that ¥ is small compared to the work-like term, 
w • A, therefore ensures that total helicity is well con¬ 
served during the redistribution of helicity by magnetic 
reconnection. 


Returning to pairing of opposite polarities of dA/dt, it 
is also straightforward to construct areas of integration 
that are smaller than the entirety of D\ and within which 
the integrated divergence term disappears. Inspecting 
Eq. (HHD, the key to this is that contours where ¥ = 0 
and curves perpendicular to Ay (components of A per¬ 
pendicular to the surface normal of D\) are generally not 
aligned. The boundary of a suitable area, S, can there¬ 
fore be built up by joining sections of T = 0 contours 
(integrand zero since T = 0) with sections of curves per¬ 
pendicular to Ay (integrand zero since A x e z • rh = 0). 
A sketch is given in Fig. O The spacing between ad¬ 
jacent ¥ = 0 contours (or points on the same contour) 
is fixed, but there is freedom to place the curves perpen¬ 
dicular to Ay arbitrarily close together. Hence, pairing 
of increases and decreases of A occurs along every curve 
perpendicular to Ay that connects two T = 0 contours. 
It is therefore seen that redistribution of helicity occurs 
in a highly organized manner between specific groups of 
field lines. 



FIG. 3. Visualization of the E 3 magnetic braid, used as the 
static part of the kinematic model with magnetic complexity, 
(left) a selection of magnetic field lines which show the braided 
nature of the field, (right) isosurfaces of horizontal magnetic 
field strength, revealing the magnetic flux rings that are su¬ 
perimposed on a uniform vertical magnetic held. Flux rings 
centered on x = 1 are right-handed and flux rings centered 
on x — — 1 are left-handed. 


V. KINEMATIC EXAMPLES 

To verify the properties identified in Section [IVJ we 
now examine a kinematic model of magnetic reconnection 
in a magnetic held with complex connectivity. The basis 
of this model is a static magnetic braid into which a time- 
dependent ring of magnetic flux is added. 

For the static component of the magnetic held, we 
use the E 3 magnetic braid detailed by Wilmot-Smith, 
Hornig, and Pontin— , which consists of a superposition 
of three left-handed and three right-handed rings of hor¬ 
izontal magnetic flux with a vertical uniform magnetic 
held. This braid is visualized in Fig. [3] and attention is 
drawn to the complex held line mapping. It was chosen 
only for convenience and any other sufficiently complex 
braid would serve just as well. 

Field line helicity is calculated subject to the gauge 
constraint Ax n \qv = A ref x n|ay, where n is the out¬ 
ward normal on dV, and 

A lei =^(-ye x + xe y ), (22) 

which corresponds to a uniform vertical reference mag¬ 
netic held, B rei = B$e z . This particular constraint en¬ 
sures that helicities are the same as for the winding gauge 
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FIG. 4. Field line helicity, A, for the E 3 magnetic braid, 
mapped on the Do surface, z = —24. Small scales are evident 
and come from the complexity of the field line mapping. 

of Prior and Yeates— , which measures winding in an un¬ 
twisted frame. It also gives H = Hr, which is zero for 
the E 3 braid. 

The map of the initial field line helicity on Do is shown 
in Fig. SJ It exhibits considerable complexity, with scales 
much shorter than those in the magnetic field (Fig. [3]) 
coming from the complexity of the field line mapping. 

Since our gauge restriction fixes the components of A 
tangent to the domain boundary and since E = 0 on 
dV when reconnection is localized within the domain, 
Eq. (|9j) implies that n x V<F|ay = 0, i.e. <£> = const on 
dV. Thus, 3>(a?o) = 3>(aq) for any field line, which makes 
the gauge terms cancel in Eq. m- Calculation of w • A 
terms is simplified too because defining w • n = 0 on Do 
and D\ gives w • A = w • A ref at field line endpoints. 
When reconnection is localized within the domain, this 
simplifies further to give 


where r is the radial distance from (x,y) = (0,0). 

Attention is drawn to the fact that provided the bound¬ 
ary condition on A is satisfied, A and dA/dt do not 
change under the allowable gauge transformations. This 
can be attributed to A obtaining a physical meaning as a 
measure of winding with reference to a fixed frame or as a 
density of relative helicity, due to the boundary condition 
on A. 

Time-dependence is imposed via a prescribed electric 
field, 

E = —Boak exp (-£- - T) e„ (24) 

which corresponds to growth of a time-dependent flux 
ring at the midplane with its center offset slightly from 
the central axis (readily confirmed by taking the curl of 
Eq. (|24D). The horizontal flux added over time changes 
magnetic connectivities between D 0 and D\, and the 


kinematic model acts as a proxy for magnetic reconnec¬ 
tion in this regard. Parameters were set as Fq = 1, 
k = 0.5, L = 2 and a = a/2/4, giving a time-dependent 
flux ring spatially smaller than the static ones, thus the 
model describes reconnection that is localized to a small 
region within the larger field. 

Compared to other models of reconnection, kinematic 
models have the limitation that the electric field is pre¬ 
scribed rather than arising self-consistently from an alge¬ 
braic Ohm’s law. The model described is nonetheless suf¬ 
ficient to verify the properties identified in Sec. IIVI which 
do not rely on any specific form for Ohm’s law. We have 
also obtained similar results using resistive MHD sim¬ 
ulations, although the kinematic model gives a clearer 
illustration because it limits reconnection in the complex 
magnetic field to a single location. 

The rate of change of field line helicity for the kine¬ 
matic model was computed as follows. First, we set 
T = 0 on Do, which also gives w = 0 there by Eq. m- 
Thus, the full evolution equation m reduces to the sim¬ 
plified form given by Eq. (USD. The values of T on D\ 
were obtained by integrating E\ \ along the magnetic field, 
in keeping with its definition in Eq. m and using the 
boundary condition T|^ 0 = 0 to fix the constant of in¬ 
tegration. Finally, the w • A term on D\ was evaluated 
using Eq. (|23|) . Note that the approach of evaluating T 
by integration of E\ \ along field lines and w on the bound¬ 
ary from VT works not only when Ohm’s law is specified 
explicitly (e.g. Eq. ©) but also when it arises implicitly 
from a prescribed electric field as in the kinematic model. 

Figure [5] shows dA/dt, and the contributions from 
the —T(xi) and w(x i) • A(x i) terms. All quantities 
are plotted mapped to Do where field lines are fixed to 
the stationary plasma by the choice T = 0 on Do, i.e. 
dA/dt = dA/dt on Do . It is immediately apparent that 
dA/dt is dominated by the w • A term, which has an 
extreme value (31.7) that is 57.6 times greater than the 
extreme value of the —T term (0.55). This finding is 
in good agreement with the analytic results of Sec. lIVBi 
which argued that the work-like term should be dominant 
given a complex magnetic field. 

The importance of complexity in the field line map¬ 
ping is demonstrated if the static braid is replaced with 
a straight uniform magnetic field (for which ^4 = 0 ev¬ 
erywhere), in which case the results shown in Fig. [6] 
are obtained. With the removal of braiding, the extreme 
values of the w • A and — T terms become 0.96 and 0.63 
respectively, giving a ratio of 1.5, i.e. neither term dom¬ 
inates strongly when the field lacks complexity. 

The kinematic example with magnetic complexity 
(Fig. [5j) also shows pairing of opposite polarities of 
dA/dt, such that field line helicity is redistributed by 
the work-like term, while the total helicity changes much 
more slowly. This too confirms our expectation from 
analytic results (Sec. IIV C[) . The analysis in Sec. II V Cl 
predicts that pairing of opposite polarities is well orga¬ 
nized. In particular, when terms are examined on D\, 
we expect to see pairing along lines perpendicular to the 
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FIG. 5. Evolution of field line helicity on Do for the kine¬ 
matic model with magnetic complexity, (top) Map of dA/dt 
(extreme value of 31.4). (middle) Map of the w(x i) • A(x i) 
term (extreme value of 31.7). (bottom) Map of the — 'f'(xi) 
term (extreme value of 0.55). Quantities are plotted on Do 
where field lines are fixed to stationary plasma; the evolution 
equation terms were evaluated at the field line’s conjugate 
point on D i, x\ — F(x,y), and mapped back to Do. 


vector potential that connect two T = 0 contours. This 
property is investigated in Figs. [7|and[8l For our choice 
of gauge, curves on D i perpendicular to A are simply 
radial lines. Moreover, the kinematic example gives a 
profile of T which is very localized on Di, T ^ 0 being 
restricted to a threadlike structure. The theory there¬ 
fore predicts that w • A (and hence dA/dt ) is organized 
in pairs of opposite polarity where radial lines cross the 
threadlike region where $ / 0. Consulting Figs. [7| and 
[HI this expectation is confirmed. 


FIG. 6. Evolution of field line helicity on Do for the kinematic 
model without magnetic complexity, in the same format as 
Fig. [5j The extreme values are: 1.34 (dA/dt, top panel), 
0.96 (w(x i) • A(x i), middle panel) and 0.63 (—T(xi), bottom 
panel). 


Note that pairing of polarities on D\ implies pairing 
on Do via the field line mapping. This is apparent when 
comparing Figs. [5] and [7[ B z = 1 on each boundary, 
and all quantities shown are defined for each field line, 
hence the figures are area-preserving rearrangements of 
one-another. However, despite the simple structure of 
pairing on D\ along radial lines, the complexity of the 
field line mapping means that the corresponding struc¬ 
ture of pairing on Do generally does not have an obvious 
analogous form. 

Returning to the topic of pairing of positive and neg¬ 
ative polarities on the top boundary, an intuitive feel 
for the exact analytic result may be gained as follows. 
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FIG. 7. Pairing of w • A polarities (hence polarities of dA/dt) 
on D\ for the kinematic model with magnetic complexity. The 
horizontal black line superimposed on the middle and bottom 
panels shows the line on which quantities are inspected in 
Fig. [8j Analytic theory predicts pairs of positive and negative 
polarities where radial lines on D\ (which are perpendicular to 
the vector potential) cross the threadlike region where T A 0. 
The prediction holds very well. 


Field lines are fixed to the fluid on Dq by our choice that 
T = 0 there, but non-ideal effects in the domain allow 
them to slip on D\. If dynamics are localized inside the 
domain, then it follows from Eq. m that motion of field 
line endpoints on D\ at any moment is along contours of 
(clockwise around maxima of and counterclockwise 
around minima of 4/). Now consider a curve on D\ per¬ 
pendicular to the vector potential, parameterized by dis¬ 
tance l along it, which connects two points where T = 0. 
In our example, the radial line considered in Fig [8] is one 



FIG. 8. Pairing of w • A on D i, along a line perpendicular to 
A ref (the section of the line y — 0 shown in Fig. [7j). The solid 
curve plots w-A (left axis) which has a bipolar character, and 
the dashed curve shows — T (right axis) which is monopolar. 
The sign of w • A changes when the direction of held line 
motion reverses. 


such line. Where d^/dl > 0, w will have one sign and 
where dT/d/ < 0 it will have the other sign. Moreover, 
the flows in the opposite directions balance because T is 
continuous and the magnitude of w is proportional to the 
gradient. The final assumption to transfer this balance 
of held line motions to the helicity equation is to as¬ 
sume that the width of the 4/ ^ 0 region is much smaller 
than the length scale over which A ref changes apprecia¬ 
bly. Then, A ref may be treated as constant to leading 
order along the curve of interest, and it is seen that the 
oppositely directed held line motions correspond to op¬ 
posite polarities of w • A (for our example, also refer to 
Eq. (j22j)). This discussion is less rigorous than the ana¬ 
lytic presentation of pairing in Sec. II V Cl but it may help 
to develop a feeling for the underlying physics. 

Finally, it is noted that the introduction of an addi¬ 
tional magnetic flux ring means that the total helicity 
of the held is changing. This happens self-consistently 
through the helicity source —2 E-B. We have already sug¬ 
gested, however, that if the preexisting magnetic held has 
a complex held line mapping, then the held line helicity 
is rearranged on a shorter timescale than that associated 
with changes to the total helicity. This is conhrmed by 
computing the flux integral of dA/dt and the flux integral 
of its unsigned value. The results are that dH/dt = 0.5, 
while the integral of B z \dA/dt\ over the domain is 11.4. 
The latter integral is 22.8 times greater than the former 
and is very strongly dominated by the w • A term. Thus, 
it is clear that reconnection in a complex magnetic held 
acts primarily to redistribute the held line helicity, and 
that this occurs much more rapidly than the total helicity 
changes. (When there is no preexisting pattern of A to 
be redistributed, which is the case for the example shown 
in Fig. [6j then the flux integrals are more similar.) 


VI. SUMMARY & CONCLUSIONS 

This paper has shown how an evolution equation can 
be derived for held line helicity, A, and examined its 
properties. For a suitable restriction of the gauge, A 
changes only if the connectivity of the held changes, ei- 
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ther via plasma motion on the domain boundary or mag¬ 
netic reconnection within the domain. Then, the evolu¬ 
tion equation becomes the sum of two terms: a work-like 
term, namely the scalar product of the generalized field 
line velocity with the vector potential at field line end 
points, and the voltage drop along the field line. 

Localized magnetic reconnection within a complex 
magnetic field was given particular consideration. This is 
relevant, for example, for magnetic relaxation of complex 
magnetic fields. Using the evolution equation for A , sim¬ 
ple scaling arguments show that evolution of magnetic 
helicity is strongly dominated by the work-like term un¬ 
der these conditions. Moreover, the work-like term has a 
property that means changes to field line helicity occur in 
paired regions of opposite polarity, which approximately 
cancel one another overall. It follows that magnetic re¬ 
connection in a complex magnetic field serves primarily 
to rearrange helicity, and has a relatively small impact 
on the total helicity. 

The idea that total helicity is approximately conserved 
on the timescale during which magnetic reconnection 
rearranges an initially complex magnetic field is long¬ 
standing. Notably, it is the conjecture underlying Taylor 
relaxations^ and has previously been justified on the 
basis that an Cauchy-Schwarz inequality implies that 
helicity decay occurs on a diffusive time scale that is 
longer than the time over which energies change, espe¬ 
cially when reconnection occurs in a small proportion of 
the total volume^. For turbulent situations, the inverse 
cascade of magnetic helicity to large scales^, where dis¬ 
sipation is inefficient, contrasts with the direct cascade 
of energy to small scales, and this has also been used to 
argue for approximate helicity conservation. This paper 
provides an alternative and complementary justification 
for helicity conservation during localized magnetic recon¬ 
nection, with the advantages that it explicitly shows the 
rapidity of helicity reorganization and provides new de¬ 
tail about its underlying mechanics. Our results also em¬ 
phasize the importance of existing magnetic complexity, 
without which the timescales of reorganization and net 
helicity change are separated at most weakly. 

Now that the evolution of field line helicity is better un¬ 
derstood, it should be possible to refine our understand¬ 
ing of relaxation via magnetic reconnection. For example, 
it is not fully accurate to say that approximate conser¬ 
vation of total helicity is the only helicity constraint (the 
origin of Taylor’s relaxation hypothesis). Rather, field 
line helicity is reorganized in a manner prescribed by the 
dominant terms in the evolution equation. The evolution 
equation may therefore reveal new constraints on mag¬ 
netic relaxation via reconnection. That possibility will 
be the subject of future work. 
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